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A UESEOI) FOR DE!DERMI]Sri£r& SHE CAUSER AND TVIST OS* A ' 
SURFACE TO SUPPORT A GITBN DISTRIBUTION OF LIFT 

By Doris Oolien. 

SUMMARY " "■ " 



A graphical method is descriljecL for finding the shape 
(caa"ber and twist) of an airfoil having an arlJitr.ary dls- 
trihu-tion of lift. The method consists in replacing ^ihe 
lifting surface and its wake with an equivalent arrangement 
of vortices and in finding the associated vertical veloci- 
ties. 

By a division of the vortex pattern into circular 
strips concentric about the downwash point instead of Tntd 
the usual rectangula:; strips, the lifting surface Is re- 
duced for each downwash point to an e^q^uivalent loaded' line 
for which the induced velocity is readily computed. ' The 
ratio of. the vertical velocity to the stream velocity is 
the slope of the surface in the free-stream direction. 

As an illustration, the shape of. the wing consistent 
with the pressure distrihution derived from the two- 
dimensional theories is found for two wings: a straight 
elliptical wing and one with 30° sweephacfc. 

Application of the method to solve the reverse prob- 
lem - finding the lift dlstr ihut ion over a given surf ace " 
is briefly discussed. 

INTRODUCTION . . - - - ; - 



Because of the effect that the. pressixre gradients 
over the 'surface of a wing have on the drag,- it would be 
of considerable advantage to be able to. specify the cam- 
ber and the twist of a wing that .would produce a desirecL" 
distribution of lift . 'Present methods' of airfoil desig'n 
depend on twor-dimensional-f low theories, which treat the 
spanwise and chordwise components of the flow independently. 
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Although a thoorotlcal troatment of threo-dlmoneloaal flow 
le givon "by Prandtl (rofcrencQ l) and calculations havo 
bocn mado for special cases (referoncos 2 and 3), no prac- 
tical proceduro for tho arhitrary lifting surface is ladl- 
catod. 

In the present paper, a method is descrilJed whereby 
the camber and the twist of a surface of arbitrary plan 
form may he determined so as to support a specified dis- 
tribution of lift. For this method, the lifting surface 
and its wake are replaced by a distribution of vortices in 
a plane. The vertical velocity induced at any point on 
the surface by the vortex system defines the slope of the 
surface at that point. Thus, the problem becomes the de- 
termination of the Induced velocities. A method is pre- 
sented for determining these veloclti«s which, by employ- 
ing chiefly graphical means, eliminates the difficult Inte- 
grations that have limited previous work. 

The substitution of a plane vortex sheet for the lift- 
ing «urfac« is analogous to the standard procedure of the 
two-dimensional thin wing-section theory (see, for example, 
reference 4, p. 87), in which- the flow about a thin, cam- 
bored section is approximated by an arrangement of vortices 
along the chord line. Inasmuch as the induced normal veloc 
Ities are assumed to be substantially tho same at the chord 
line and at the airfoil, the ratio of those velocities to 
the fres-stream velocity gives the slope of the camber line 
In three-dimensional flow a reference plane is assumed, so 
situated that the airfoil may be considered to be a slight 
deviation from it. Upon this plane the plan form and the 
pressure distribution of the surface are projected. As in 
two-dimensional treatments, the slope is calculated in the 
froe-stveam direction. 

The substitution of a vortex sheet for a lifting sur-." 
face is discussed at some length by von Ki^rmAn (reference 
5, p. 15). In the application of the method, however, the 
vortices are generally assumed to have a rectilinear dis- 
tribution. Even with this limitation, the evaluation of 
the integrals involved in finding the induced velocities 
presents considerable difficulty. (See oh. IV, sec. 15, 
of reference 5, where the formulas are dbveloped for a 
rectangular wing.) The integration is greatly simpliflod 
by the introduction of polar coordinates so chosen that 
the elemonts of integration are circular strips concontrio 
about the point at which tho downwash is to bo f^ound. Ap- 
plication of this mothod is not restricted to any partlcu- 
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lar ar raagemen-b of vortices or to any specific form of the 
surface . ' . . - 

Although the procedure deecrihed is for the detormina' 
tion of the surface that will fit a req.uired pressure dis- 
tri'butioa, it may also "be adapted to effect the reverse 
analysis, that is, to find the pressure distribution over 
an arlJitrary surface. When the surfaces that fit several 
assumed pressure distrihutions have "been found, the re- 
quired surface can be built up by a linear combination of 
th.ese solutions* If a reasonably close distribution can 
be assumed as a first approximation, finding the surface 
to fit the assumed distribution will indicate the' manner 
in which that distribution must be corrected to fit the 
given surface. 

DBTBEMI STATION OF TH3 VOHa?EX PATSEEH 

FBOM THE FBESSUHE OISTHIBUT lOK 



The vortex pattern is obtained by integrating the 
chordwl.se pressure distribution back from the leading edge 
at several stations along the. span. The cipculatioh or 
vortex strength r will be shown to be proportional to 
this integral; the lines connecting the points where the 
values of the integral are equal therefore^ define the vor- 
tex lines. 

In figure 1, a distribution of lift is arbitrarily 
specified for a tapered wing in straight flight. The cor- 
responding vortex lines are drawn on the plan form of ^he 
ving and in the wake in figure 2 to show a typical vor"^ex 
pattern. 

The demonstration of the relation between the pres- 
sure distribution dnd the vortex pattern is given in the 
following paragra.phs. 

In the replacement of a lifting surface by a vortex 
sheet, the assumption is made that the pressure increments 
due to the presence of the airfoil in the stream are equal 
and opposite on the upper and lower surfaces, as Vould be 
true in the case of a thin plate at a small angle of at- 
tac-k. The substitution is still admissible in the calcu- 
-latioa of lift when the thickness. is not negligible, be- 
cause the difference between tie velocities on the upper 
and lower surfaces at any position and not' the magnitude 
of each increment determines the lift at that point. 
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Thus, let AV^ and AT^ represent for any point the 

local velocity increments on the upper and the lower sur- 
faces due to the plate or to the rortices that are eq[ulva- 
lent to it; let u^ and uj be the magnitude of their 

cempozients in the direction of the free-stream velocity T 
let v^ and Vj he the components normal to V In the 
plane of the wing; and let w^^ and w; "be the component s 

normal to the wing. Then, the' pressure on the surfaces 

would "bo ........ 

Pu - I 1 + 

= i p (V^ + 2Tu^ + u^^ + v^* + w^^) 

= i p (7^ + 2VUu) 

no^le c ti-ng second- order effects ; nnd, similarly, 

Pl = I P (V^ + 3Vuj) 



The resulting lift per unit area is then the difference in 
pressure, or 

= ^ P2T (u^i - Ml) (1) 

The derivation of the equivalent vortex pattern fol- 
lows directly from equation (l). If dgT represents the 
element of circulation around a small length ds, paral- 
lel to Xi over which the velocities u^ and ii^ may 
"be considered constant, the following relation holds: 

d^r = (u^ - uj) ds (3) 
from which equation (l) may be rewritten 

Ap = pV |i: C3) 
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Thus, the lift at evory Doiut is proportional to »' or 

the cross-stream component of the vorticity. Consider now 
a narrow strip of varying width just "behind the leading 
edge of the vortex sheet, suoh that / Ap ds is constant 

along the strip. Such a strip would represent a vortex 
element of strength / Ap ds. A second vortex element 
could he defined in the same way to lie just hehind the 
first. ffrom eq_uation (3), however, 

Equation (4) defines the function P for any point on 
the lifting surface as the total circulation ahead of the 
point.* Thus, the "boundarios of the vortex element are 
linos of equal F. 

If the vortex elements are reduced in size and in- 
creased in number to form a continuous vortex sheet, their 
pattern is Indicated, as in figure 1, hy the contour lines 
of the function P. In order to satisfy the Eutta 6*ondi- 
tion that there he no pressure difference across the trail- 
ing edge, these lines must leave the wing paraliel't"© the 
stream velocity, as shown by equation (3). They then fol- 
low the streamlines down the wake. The drawing of these 
contour lines from the integral of the pressure distribu- 
tion is the first step of tJie procedure for finding the 
camber and the twist of the lifting surface', "" ^ 



DETERMIHATIOIT OF THE INDUOED Y3ETICAL VELOCITIES 



The induced vertical velocity at any point of the 
plane is obtained by integrating the Biot-Savart oquation 

*Th6 function T is numorically equal to the differonoo 
in velocity potential between the surfaces, for, at any 
point, u^ = o0^/ds and u; = S0i/Ss, where 0^ and 0; 

aro the potential functions over the upper and lower sur- 
faces of the airfoil, respectively . From equation (2), 

ds OS 3s 

or 
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over the entire vortiex pattern Just de scribed. The result- 
ing double integration is reduced to a single integral by • 
a relatively simple graphical procedure if a Byatem of 
polar coordinates having its origin at the point is em- 
ployed. *■ 

Let it be required to find the downwash w at a point 
P of the vortex pattern just described. Consider a small 
sector of the plane included between two radii from P, 
the angle between the radii being d\J/ (fig. 3). At any 
distance r from P, the radii will cut off a small 
length dl of a certain vortex element. If the width of 
the element in the radial direction is dr, the strength 

of the vorte-x element is - ^ dr. Then, by Biot-Savart ' s 

dr 

rule, the downwash at P due to the small length of vor- 
tex will be 

dw = - 1^ dr dl sin ^ (B) 

%TT r or 



where p is the angle between the vortex element and the 
radius. The length 6.1 sin 3 is the projection of dl 
on the circumference of the circle of radius r around P. 
so that 

dl 'sin P = rd\|/ 

and 

dw = - -r^ ^ dr d\lf (6) 
4TTr dr 

The total downwash at P would be obtained by integrating 
(6) over 0 ^ r >» and through. 360 of y\r . Thus, 

CO SIT 



In order to evaluate the double integral of eq.uation 
(7), the relation (Liebniz's rule) 

3TT STT 



f If "i^l' = y r dM; (8) 
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is used to reduce the first integration to a graphical 
procedure, 

The first step in evaluating v is to draw on the 
vortex pattern circles of various radii aTsout P. Around 
any one circle (see fig, 4) the function P will tafco on 
values indicated "by the inter soot ions of the circle with 
the contour linos of F . If those values of F are .plot- 
tod against the angle ^ -moasuredt let us say, from 
the free-stream direction - gi?aphical integration will 
Sir 

give / r d>I^ for each circle. It is somowhat more coa- 
o 

venient to plot F against '^/2rn. Then the integral will 
"be the average value of P around the circle, designated 
in the usual way "by 

(9) 

Then F is a function of r. (See fig, 5.) I'rom ecLuation 
(8), 

air 




J 8r 



so that 



d\|f = Sit 1^ (lO) 



(11) 




I'or the evaluation of equation (ll), first plot F 
against r, as In figure 5. The curve will approach an 

asymptote, which is readily found from the ezpressior. for 
the load curve across the wake. Thus, when r is large, 

= and eq^uation (9) reduces to . 

■b/2 



■if 



r(y) dy (12) 
-■b/2 ._ 

where h/s is the semi span of the wing. If the area un- 
der the load curve across the wake is A» then for very 
large values of r ' . 

r = -4- (13) 

2tTr 
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The plot of r(r) should he carried out to a value of r 
such that the curve approaches this asymptote within the 
accuracy of the work. 

The load curve of figure 5 is a typical one. The fol 
lowing method has hoen found particularly suited to the 
determination of the downwash at the origin of such a 
curve. The first section of the curve (designated lay l), 
starting with zero slope, is approximated up to the infloc 
tion point, or to a point *q somewhat ahead of it, hy 
an expression of the form T' = a^ - eigr^ - a^r* . Addition 
al tiprms might he used hut are generally not necessary. 
Tho downwash due to a curve of the general form 
r = ao - aiir°-, (n > l) , at r = 0 is given hy 

w = - ^ ^ ^ 2. (^B^"^ - ^a'^"^)' where r^ and are t; 

end values of r for .the interval over which the curve ex 
tends* The downwash contrihuted hy the first section of 
the load T = ftg - a^r® - a^r* is therefore 

w(l) = - fsaaro + ^ a^ro^) (14) 

The part of the curve iiBHiediately following r = 
has a critical effect on the value of the downwash, at the 
same time being . usually too irregular to bo approximated 
for any distance by a simple algebraic expression. It is 
therefore advisable to proceed as in numerical integration 
dividing the curve. into a finite number of small sections 
and' considering each section of the curv« to have a simple 
mathematical expression. Because the effectiveness of the 
variation of load depends on its closeness to the downwash 
point, the intervals are taken in geometric rather than 
arithmetic progression. Thus, the abscissas are r^, krQ, 
k^rQ, . . , etc, where the ratio k is a number, usual- 
ly between 1 and 2, determined by the sisse of the Inter- 
vals required for accurate representation of the curve 
immediately following ttq. The usual procedure now would 
be to assume the curve to ba a straight line over each 
small interval but when the curvature is largely in one 
direction, as it is in these curves, this assumption in- 
troduces a small but cumulative error, which may amount 
to 10 percent or more in the total. The following method, 
which fits the curve with a succession of parabolas, is 
found to give very good accuracy with no increase in com- 
putation. The method is best--presented in tabular form: 
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an 



311+1 



an 



(Ooarpute) 



AT 
an 



S-1 



'•a 



^2N-1 



r -r 



An 



AH 



aT-2 



22r-2 



An 



n=o n=o 



231 + 1 



log k 

^1 = ^ — kijr 



10 



NACA Technical Note ITo. 855 



In the following ta^ble ^.re given the constants Zq, , 

and — 7— r for several convenient values of fc. Since 

k(k-i; 

the choice of k is not critical, t"h« values included 
should serve without interpolation. 





Ko 


Ki 


1 

k(k - 1) 


1.02 


0.01000 


0.01000 


49.02 


1.05 


.02582 


.02420 


19.04 


1.07 


.03667 


.03340 


13 .35 


1.10 


.05325 


.04690 


9^091 


1,15 


.08222 


.06827 


5.797 


1.30 


.1127 


.08840 


4.167 


1,25 


.1446- 


.1074 


3.200 


1.30 


.1780 


.1255 


2.564 


1.50 


.3246 


.1891 


1.333 



The value of the ratio k that- will give sufficiently 
small incroments of r where the slope is large will proh- 
ahly Tje found to "be smaller than necessary after the ourve 
has become more regular and the distance from the origin 
larger. The computations may he interrupted here and the 
downwash w(ll) contrihuted hy the second section may he 
cnlculated. Then, using the last abscissa rg^ as a new 

starting point and a larger value of k, compute in the 
same way the downwash w(IIl) due to the remainder of the 
curve to a point where th« difference between the curve 
B.nd its asymptote is negligible. 

ffor the portion of the curve extending to infinity 
(section IV of- fig. 5), the previously determined 'asymp- 
tote is used and the downwash found analytically. Thus, 
from a large value E or r out to infinity, 

oo 
R 

from equation (lO), and if dT/dr is found from equation 
(13). 
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The downwash at Pis then the sum 
w(I) + w(ll) + w(III) + v(lV). 

It is interesting to note that the three-dimensional 
problem has "been reduced to one of two-dimensional flow,' 
as mar he seen "bjr replacing F with its original expres- 
sion (9) in eq.uatlon (ll), which may then be written 



00 



d / r d\if 
^ — dr 



w(o) = / (16) 

4Tr (r - 0) 
o 

Equation (16) is recognized as the ordinary formula for 
the in'duced normal velocity with the load expressed in the 
form of a definite integral and suggests that the ^irst 
integration (e_xcept that the factor 1/2tt was introduced) 
was equivalent to concentrating all the vorticity around 
each circle at a single point at the distance r Along 
a line of infinite extension from' P. The loaded line of 
figure 5 may be considered, except for the factor . l/Zn, 
to be the eq.uiyaleat of the original lifting surface. 
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The method will be applied to check the elliptical 
distribution of lift cohventionally assumed for an uncam- 
bered elliptical wing. Shis distribution, arrived at by 
combining the two-dimensional theories, does not take ac- 
count of sweepback or stagger of the lifting elemements. 
In the applicati.on of the present method, two cases of ol- 
liptloal chord distribution will bo inyestigatqid, one with 
a straight SO-per ce nt -chord line, and'^no^ with tho 50- 
percont-chord line swept back about 30 , the sections re- 
maining parallel to the plane of syiSmd'try.' The calcula- 
tion of the vertical velocities over the swept-back wiiig 
will bo carried out" herein in some detail in order to il- 
lustrate tho method. The downwash will be found at sev- 
eral points along the three-quarter-chord line. Tho down- 
wash at the three-quart er— chord point of a "section is of 
particular intorost because, from tho thin wing section 
theory (reference 6, p. 82), the effective angle of attack 
of the section is given by the slope of tho camber line 
in the ne-ighborhoad of that point, if the camber is approxi- 
mately circular. ■ ~' 
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For the chordvise lift diBtrnJutlon, the two-dlmonBi on- 
fil flow around a flat plate as givon "by the thin wing sec- 
tion theory is assumed. This lift distrihution and the 
circjilation function obtained by iategrnting it aro shown 
in figure 6. In this oase the mathematical expreesion for 
the plresBure difference is known and P can "be found ana- 
lytically, ■ For convenience, the uni-bs have been so choson 
that the maximum' value of P is 1.0. If, further, all 
lengths arc ezprossed in terms of the somispan (measured 
perpendicular to the plane of symmetry) as the unit of 
length, t he tota l circulation (T at tho trailing odgo) is 

simply yi - y®, where y is tho distance spanwiso from 

tho center line. The resulting contour lines of P ap- 
pear, for the swept-back wing, as in figure. 7, The points 
at which the downwash . will "bo found are also shown. 

In order to find tho downwash at point B, clrclos 

spaced as shown in figure 8 are drawn about this point 
and tho Values of P indicated by the contour linos inter- 
ooptod by each ..circle are plotted against tho angular lo- 
cation ^ of the points of inter sixct ion. Of the ourvos 
for r(\|/; corresponding to the circlp;s of figure 8, fivo 
typical ones are shown in figure 9. Tho curves dosignatod 
for r = 0,067 and r = 0.20 are characteristic of cir- 
cles, close to tho downwash point; the curve for r = 0.216 

includes a point (^/2v = 0.883) at which tho circle is 
tsngent to tho leading edgo of the wing; tho circle of 
radius 0.733 Ilea partly ahead of the wing where P = 0, 
partly in the wing, and passes through the wake; the cir- 
cle of radius 1,60 traverses tho vortex pattern across 
tho wake only. 

The results of integrating these curves are plotted 
against r in figure 10. Since the cinculation acrosB 
the wake la ^ 1 - y^ , the integral iT i-s tt/2 and the 
asymptote for the curve of Vir) is, from equation (l3), 
r = l/4r. The down\irash is now calculated as follows: 

It is found that the polynomial 

r = 0.8335 - 2.70r^ - 6.57r* 

fits the curve of T through r = 0.20. Then the down- 
wash contributed by the section from 0 to 0.2 is, from 
equation (l4) 

w(l) = 0.575 



KAOA Technical ITote No. 855 



13 



Section. II, with, k small, is taken to include Isoth 
inflection points (there is a third inflection point at 
1.78, hut its effect is negligihle) . Tahle I shows the 
calculations for this section and for the succeeding one. 
Section III extends to r = 2.433, wher6 'T = 0,108 and 
l/4r = 0,103. The downwash induced hy the asymgtoto from 
r =3.433 out to infinity is only 1/33(2.433) =_0.011, 

ahout 1 percent of the total due to the curve of F, and 
the S^percont error in the ordinate may therefore ho nog- 
lectod. The total downwash at -B is then 

w(I) + w(Il) + w(III) + w(IV) » 1.156 

The downwash at points A, C, and D iB_found in the 
same way, I'igure 11 shows the curves of T for those 
points. The asymptote P = l/4r is, of course, common to 
all the curves for this wing, 

The downwash is plotted against the spanwisp location 
of the points in figure IS,, Tho quantities Tjag^x 
h/2, heretofore asstimed to he unity, are included to make 
the result nondimensi onal , This curve of the downwash or 
vortical velocity at the three-quarter-chord lino, sinco th 
slope of the surface is w/V, ^ is a measure of the amount 
of twist the wing must have to sustain the. assumed distri- 
hution of lift. If the wing were actually flat, which was 
the premise in deriving the distrihution from the two- 
dimensional theories, w would he equal to 1.0 all over 
the surface, so that the deviation of the curve from the 
line w = 1 indicates the amount hy which tho two-dimen- 
sional theories are in error when applied to three- 
dimensional flow. The discontinuity in vorticity at tho 
center line gives rise to a discontinuity in the downwash, 
which goes to infinity everywhere along the center sec- 
tion. This result indicates that the assumed condition, 
in which the vortex lines hend to form an angle, cannot 
exist in practice. The vortex lines actually would he 
rounded off and the load at the center reduced "below that 
of the adjoining sections. 

Tho corresponding curve for the straight elliptical 
wing is shown in figure 13. In this case, the downwash 
Was computed also at points along the quarter-chord line. 
The slope of the surface was found to he less than the 
slope at tho corresponding three-quarter-chord points. 
The positive oamher thus indicated is very small at the 
center hut increases sharply near the tips. 
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It is expooted that the method outlined herein will 
bo ospocially useful in investigating the nature of the 
flow near the wing tip, where two-dimenai onal ajpproxlraa- 
tions no longer can "be applied. The more accurate threo- 
dimenelonal treatments available are also unsuitable for 
this purpose becauso the calculations fail to converge at 
the tips. The method of calculation described In this pa- 
per presents no particular difficulty in these regions. 
As a test of the accuracy obtainable, the induced down- 
wash was computed at a point near the edge of a circular 
plate in nonlifting: potential flow and was found to check 
almost exactly with the known solution. 

Kevertheless , results obtained for the tips have 
qualitative rather than quantitative value (except for low 
angles of attack), The validity of the theory is actually 
limited by the existence of strong tip vortices, which may 
cause the vortex sheet to curl up out of the plane in 
which it is assumed. to lie. On the other hand, the high 
concentration of vortlcity associat-ed with this effect 
adds appreciably to the drag of the wing, so that even a 
general indication of such a concentration of load at the 
tips is of-^alue. It should be possible to design, by 
the use of the present method, a wing tip that would avoid 
this effect by providing a fairly gradual tapering off of 
tho load spanwise. At the same time, a favorable chord- 
wise gradient could be- specified. Thus, it appears likely 
that an optimum tip for low drag could be deduced. 
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TABLE 1 
DOWSWASH COMPUTATIOITS 



r 


r 


4r 


^an 




Section II: k = 1.07 


0.200 
.214 
.229 

.245 
.26 2 

.280 
.300 

.321 
.344 

.368 
.393 


0.715 
.690 
.657 
.636 
.620 
.606 
.593 
.582 
.569 
.560 
.550 


-0.025 

-.033 
-.021 

-.016 
-.014 

-.013 
-.011 

-.013 
-.009 

-.010 


-0.125 
-.092 
-.053 
- . 033 
-.026 


-0.165 
-.070 
-.050 
-.043 
-.029 




Totals: 


w(ll) 0.330 


-.329 


-.357 




Section III: k s 1,20 ' 




0.393 
.472 
.566 

.679 
.815 
.978 
1.173 

1,408 
1.790 

2.028 
2.433 


0.550 
.524 
.497 
.453 
.379 
.317 
.268 
.219 
.153 
.131 
.108 


-0.026 

-.027 
-.044 

-.074 
-.062 

-.049" 
-.048 

-.06& 
-.022 

-.023 


-0.066 
-.078 
-.076 
-.042 
-.012 


-0.069 
-.131 
-.060 
-.056 
-,013 




Totals: 


w(lll) = 0,250 


-.274 


-.329 

1 




Figure 1.- An axbltxazy dlstzlbutlon of lift assumed for a tapered wing In 
straight flight, from whloh the vortex lines of figure a were 

derived. 




Figure 2. 



_ Contour lines of olroulatlon funotion, or vortex pattern, obtained by integrating 
the pressure dlstzlbutlon of figure 1. 
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Figure 12.- Vertical velocity w at the three-quarter chord line of the swept- 
back wing with distribution of load calculated by two-dimensional 
flow theories. 
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Fig. 13 




Figure 13,- Vertical velocity w at points on stralBht elliptical wing, 

A K 6, with distribution of load calculated by two- 
diaenslonal-f low theories . 



